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Abstract. Recently an additional technique was applied to investigate the properties of kHz QPOS, i.e. the analysis of the 
distribution of frequency ratios or frequencies themsel ves. I review the res ults of such work on the data from Sco X-1: 
lAbramowicz et all i2003h . which was later criticized bv JBelloni etaill2005 h. I find that the findings of the latter paper are 
consistent with results presented earlier: kHz QPOs cluster around the value corresponding to the frequency ratio of 2/3. I 
also discuss the random walk model of kHz QPOs and possible future observations needed to verify it. 
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1. Introduction 

The nature of kHz quasi periodic oscillations continues to 
puzzle astrophysicists. In order to unveil the nature of these 
phenomena different quatities are considered. For example 
the difference of the two kHz QPO, which in the beat fre- 
quency mo del is constant was analyzed . The disk resonance 
model e.g. Abramowicz and Kluzniak (2001) predicts that 
the ratio of the two kHz QPO frequencies should assume a 
set of discete values. In order to test it we analyzed the statis- 
tical p roperties of KHz QPOs in Sco X-1 ('Ab ramowicz et all 
l2003h . We have found that the frequency ratio of the lower to 
higher QPO cluster around 0.688, a value close to the canon- 
ical 2/3 seen in the QPOs from black h ole candid ates. This 
discovery has recently been disputed by Belloni et al. (2005) . 

In this paper I will discuss the similarities and differences 
of the two analyses. In section 2 I present an overview of 
different methods used to compare distributions and in sec- 
tion 3 I summarize the results of Abramowicz et al. (2003). 
I compare our results with the work of dBelloni et al.L 120051) 
in section 4. Section 5 contains summary and a discussion of 
the possible future observational tests. 

2. Comparison of distributions 

Inferring the shape of an underlying distribution given a set of 
observed quantities is a diffcult task. In a number of cases we 
have no prior knowledge of what the shape of the distribution 
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should be. Very often we do not know if the measured quan- 
tity is fundamental from the physical point of view. However, 
one can start by making simple assumptions about the data 
and testing them. 

There is a number of ways to compare the shape of two 
distributions e.g. a model and an observed one. The chi^ test 
comparing a binned observed distribution with a model are 
not very useful. Such comparison suffer from a number of ar- 
bitrary assumptions about the widths of the bins which makes 
the results not reliable. A different set of statistical tests that 
do not require binning of the data is based on comparing cu- 
mulative distributions. One can introduce various metrics to 
measures distance between two cumulative distributions and 
calculate the probability distribution of such distance. A very 
useful test using this approach is the Kolmogorov-Smirnov 
test in which the metric is defined as 

i?K5 = max|(Ci -C2)|, (1) 

where Ci are the cumulative distributions. The sensitivity of 
such test is limited, yet it serves as a convenient tool for com- 
paring distributions and assessing goodness of fit. There exist 
modified and more sensitive versions of this test. For example 
the Fiszer-von Mises tests uses the metric defines as 

DF^.z = J (Cl -C2)^ (2) 

and is therefore much more sensitive to differences in overall 
shape of two distributions. However calculating the proba- 
bilty distribution of D pisz is a much more demanding task/ 

Any such test can be used for parameter estimation. Sup- 
pose that the model distribution Cm {Pk ) is a function of sev- 
eral parameters pk- We can then maximize the probability of 
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Fig. 2. The frequency ratios used by lAbramowicz et all 
("2003) are represented by the vertical lines. The solid line 
corresponds to the best fit single Lorentzian. 



Fig. 1. The binned distributions o f the lower and upp er QPO 
frequencies from the data set o fvan der KHs et all {^97). 
Note that the two distributions are similar but not identical. 



the model given the observational set of data. In this case 
however the most sensitive method to use is the maximum 
likelihood. In this approach we can calculate the probability 
of the model M given the data X : (xi): 

C{M\X) ocH^ixi). (3) 
ax 

This likelihood function can be maximized with respect to the 
model parameters to yield the best fit model. The goodness 
of fit of such model can evaluated using e.g. the Kolmogorov 
Smirnov test. In order to find the probability distribution of 
the parameters one ususally assumes that it is proportional to 
the likelihood function. 



3. Summary of the Sco X-1 analysis 

In our paper we have addressed the question whether there 
is any evidence of clustering of the frequency ratios of kHz 
Q POs in Sco X-1. To do tha t we used a set of data published 
bv lvan der Klis et aD ( ll997l) . This set contained 39 frequency 
pairs and was used by the authors to prove that the differ- 
ence between the two frequencies is not constant. We have 
first verified the hypothesis that the distribution of frequency 
ratios is uniform using the Kolmogorov-Smirnov test. This 
probability is 

-Puniform = 0.00047, (4) 

which led us to consider the question if the frequency ratios 
cluster We have no a priori knowledge of the shape of the 
distribution we are looking for as there is no physical model 
that we could refer to. Therefore we decided to model the 
shape of the distribution with a Lorentzian profile normalized 
on the interval (0,1). We have used the maximum likelihood 
method to find the best fit which was a Lorentzian centered 



on To = 0.688 with a width of A = 0.0151. The Kolmogorov- 
Smirnov test provided a confirmation that the data could have 
been drawn from such distribution, since the probability was 



Lorentzian 



= 0.678 



(5) 



We have also shown that there possibly is another peak in the 
data at r w 0.78, yet is was not significant. 



4. Comparison with the later study 



In a recent paper lBelloni e t al' ^2005^ criticize the results ob- 
tained above. There is a relation between the two kHz QPO 
frequencies, and therefore it is equivalent to analyze the dis- 
tribution of one or of another frequency. The relation between 
the two kHz QPO frequencies is cle a rly vi sibly and could 
already be noted in Ivan der Klis et al] ( Il997l) . The statement 
would have been true if the relation was exact, however this is 
just a correlation and not a strict one to one correspondence. 
One can verify that by checking the consistency of the distri- 
butions of the lower and upper frequency with the hypothe- 
sis that they are distributed uniformly, in an interval between 
the lowest the highest frequency in a gi ven set. For the uppe r 
frequency we obtain the same result as i Belloni et al.' feOOSi). 
i.e. Puppor = 0.017, while the analysis of the distribution 
of the lo wer frequency gives -Piowcr — 0.0012. In their con- 
clusions 'Bellon i et alJ (12005') write that the "distribution of 
QPO frequencies used by Abramowicz etal 2003 differs from 
constant can only be rejected at 2.4(t level" This statement 
is only true for the distribution of the upper frequency and 
not for the lower. Moreover we ne ver considered the d i stribu- 
tions of individual frequencies in lAbramowicz et al] (l2003l) 
but their ratio. The distribution of the frequency ratios differs 
from uniform at a level of 3.5(7. Therefore the statement that 
"the results of Abramowicz et al. 2003 are not statistically 
founded" is false. 

iBelloni et alJ ( l2005l) analyze the distribution of the up- 
per kHz QPO in Sco X- 1 using a larger data set. They found 
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that the distribution of this QPO is not compatible with a uni- 
form level and the fit the distribution using three gaussians. 
It is interesting to note that in all cases there is a peak in 
the distributions they obtain corresponding to the 3/2 ratio 
of the frequencies assuming the linear relation between the 
two frequencies. The distribution of the frequency ratios can 
be transformed to the distribution of an individual frequency 
using the lines relation. Using the propagation of errors we 
estimated that there should be a peak in the distribution of 
upper QPO in Sco X-1 at the frequency 
t^peak = 940 ± 25Hz. (6) 
iBelloni et all ilOO^ found a peak at 432.5 ± 1.5 Hz in com- 
plete agreement with our results. It is extremely interesting to 
find that the distributions of QPO frequencies in other sources 
analyzed bv Belloni et al.. C200 5) also show clustering and 
peaks at the location consistent with the value of 3/2, when 
the linear relation between the two kHz QPOs is taken into 

acc ount. 

iBelloni et alJ ( l2005l) propose an explanation for the clus- 
tering of frequencies as a simple randon walk. In their model 
they start from 700Hz and let the frequency wander in fre- 
quency in steps of ±6 Hz. They show that this model quali- 
tatively explains clustering of frequencies. This model is not 
defined in detail, as it has not been clearly stated for how 
many steps the frequency is allowed to wander, or if there are 
any boundaries outside of which the frequency can not drift. 
Moreover in a real observation the starting point depends on 
a particular moment of the observati on. The distributions of 
th e upper kHz QPO of Belloni et al. (2005) is similar to the 
of Ivan der Klis et all |,1997), and the Kolmogorov-Smirnov 
test shows that the probability that they come from the same 
underlying distribution is 22%. One can ask a question how 
probable is it to obtain two such similar distributions in two 
realizations of a random walk. We have preformed a simu- 
lation and found that the probability of obtaining two such 
similar distributions varies between 10^'' and 10^^ depend- 
ing on the specific prescription of the random walk. 

5. Summary and future observations 

We have shown that the analysis of the d istribution of t he up- 
per kHz QPO frequencies in Sco X-1 bv'Bell oni et"al] ( l2005l) 
leads to the results in to tal agreement with the results of 
lAbramowicz et alJ ll2003l) . The pe aks in the distributio n of 
upper kHz QPO frequency found bv'Bell oni et alJ (l2005l) cor- 
respond exactly to the es timates based the earlier results of 
lAbramowicz et 'ai] ll2003l) . The deviation from a uniform dis- 
tribution of the frequency ratios is at the 3.5(7 level, and not 
2.4(7, which corresponds t o the dis tribution of the upper kHZ 
QPO as claimed by Bell oni et alJ (^05). the distribution of 
the lower kHZ QPO differs from uniform at a 3.2(7 level. 
The difference comes from the fact that the linear relation be- 
tween the two frequencies is not exact but has s ome spread. 

The random walk model proposed by 'Belloni et alj 
(I2OO5..) reproduces qualitatively the distribution of upper KHz 
QPO. The problem that it faces comes from the fact that it is 
very difficult to reproduce a similar distribution in two sep- 
arate observations, while the set of data in van der Klis et all 



1I1997I) is consistent with the set of lBelloni et"ai] (|2005). This 
model can be verified with a long term monitoring of sources 
with kHz QPOs and a uniform sample of their behavior If the 
frequencies, and frequency ratios really do cluster on some 
values such behavior should be reproducible in consecutive 
long observations. In the random walk model the peaks in the 
distributions of frequency ratios will wander from one obser- 
vation to another Such study must be treated with caution as 
we have no information on the fundamental scale of variabil- 
ity of kHz QPOs. 

Furthermore in my opinion there is no preferred way of 
analyzing the data. There is no reason a priori to say that the 
frequencies or frequency ratios, their sum or difference is in 
any way a better quantity to analyze until we have a phys- 
ical model to consider. Abramowicz et al. ( 2003) based the 
analysis on such a physical model in which the frequency ra- 
tio plays a fundamental role. The analysis bv IBelloni et alJ 
( 2005 ) is phenomenological and does not refer to any under- 
lying physical model. 

Finally there is the issue of the con sistency of the data 
with the resonace theory ( Abramowicz and Kluz niak. 200 li). 
The theory predicts a preference for the 2/3 ratio of the fre- 
quencies but is not saying anything yet about the shape of the 
expected distribution around this value. Such shape should 
probably depend on the details of the interaction in the inner 
disk. At this point we can only say that the value of the ratio 
r = 2/3 lies within the width of the assumed shape of the 
distribution of frequency ratios in Sco X-1. This observation 
provides an additional hint to consider the resnance model in 
more detail. 
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